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Magnetic skyrmion, a swirling spin texture, in chiral magnets is characterized by (i) nano-scale size (~lnm 
- lOOnm), (ii) topological stability, and (iii) gyro-dynamics. These features are shown to be advantageous 
for (a) high-density data-storage, (b) nonvolatile memory, and (c) ultra-low current and energy cost ma¬ 
nipulation, respectively. By the numerical simulations of Landau-Lifshitz-Gilbert equation, the elementary 
functions of skyrmions are demonstrated aiming at the design principles of skyrmionic memory devices. 


1. Introduction 

High density, high speed, low-energy cost nonvolatile memory devices are regarded as 
an essential element of the next-generation electronics. Many candidates have been 
explored up to now in this respect, such as Magnetoresistive Random Access Memory 
(MRAM), Phase change Random Access Memory (PRAM), Resistance Random Access 
Memory (ReRAM), etc. 1-6 In addition to these existing memory devices, we propose 
in this paper the memory function utilizing skyrmions (Sks) 7 as the new promising 
candidate which has many superior characteristics. Figure 1 shows schematically the 
magnetic configuration of an Sk in magnet. It is a vortex-like structure in the ferro¬ 
magnetic (F) background with the magnetic moments pointing down at the core and 
swirling in the intermediate region. Because the winding plane is perpendicular to the 
radial direction in Fig. 1(a), it is called Bloch Sk. It is noted that an Sk has a finite size £ 
typically of the order of l-100nm, in contrast to a vortex or a meron, 8,9 which is “half’ 
of an Sk with the swirling in-plane components of the magnetic moments extending to 
infinity. This size £ is determined by the ratio of the relativistic Dzyaloshinskii-Moriya 
(DM) interaction 10,11 D and the exchange coupling J as £ ^ ( J/D)a with a being the 
lattice constant. The DM interaction is allowed in non-centrosymmetric chiral magnets 
such as MnSi, 12-17 Fei-^Co^Si, 18-23 FeGe, 24-27 and Mni^Fe^Ge 28 and Cu 2 0Se0 3 . 29 In¬ 
tuitively speaking, DM interaction prefers the twisted spin configurations while the 
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Fig. 1. (Color online) Schematic illustration of a skyrmion. (a) Bloch skyrmion. Prom the center to 
the rim, the rotating magnetic moments are always perpendicular to the radial vector, (b) Neel 
skyrmion. The rotating plane of magnetic moments is perpendicular to that of Bloch skyrmion. The 
color map (right) specifies the in-plane components of the magnetic moment and the brightness of 
the color represents the out-of plane component. 


external magnetic field stabilizes the F one. To reconcile these two interactions, the 
spiral magnetic configurations are preferred and as the special form among them, the 
skyrmion crystal (SkX) is realized in the intermediate region of the external magnetic 
field. In the multilayered magnet, the uniform DM vector is anticipated to exist along 
the interface of the layers, and hence a Neel Sk shows up 30 which is schematically shown 
in Fig. 1(b). These Neel and Bloch Sks are in the same topology class as discussed below. 

The topology of Sk is characterized by an integer N s k called skyrmion number which 
counts how many times the direction of the magnetic moments n r wrap the unit sphere. 
Both Bloch and Neel Sks have the same N s k- This topological integer does not change 
as long as n r remains slowly varying compared with the lattice constant a, i.e., n r is 
viewed as a continuous function of the spatial coordinates r. This provides the stability 
of an Sk as an emergent particle composed of many magnetic moments. For example, 
an Sk created in the F background has a very long lifetime even though it is a meta¬ 
stable state. In other words, the discontinuity of the magnetic configuration must be 
introduced to create or annihilate Sks, which we call “topological surgery”. 

As for the dynamical properties of the skyrmionic systems, it is noted that the 
effective magnetic field is associated with the Sk. Namely, the solid angle subtended by 
the spins leads to the Aharanov-Bohm effect and acts as an emergent magnetic field 
both for the Sk itself and conduction electrons coupled to the Sk. 7 The equation of 
motion under the spin polarized electron current j s (of which sign is opposite to the 
electric current j) for the center-of-mass motion of an Sk reads 

+ G x (j s - v d ) + k ( av d - / 3j s ) = - W. (1) 

Here (A, Y) is the centre of mass coordinate and v,i = (A". Y) is its velocity. The mass 
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m of the Sk originates from the deformation of the moving Sk compared with the 
static solution. Recently, this mass m has been studied extensively, and it is found that 
m is almost zero when the motion is driven by the electric current. 31 The gyrovector 
G = (4irN s k)e z (e z : unit vector along ^-direction) describes the gyro-dynamics of the Sk 
analogous to that of a charged particle under magnetic field. A dimensionless constant 
k is of the order of unity, a denotes the Gilbert damping constant, 3 represents the 
coefficient of the non-adiabatic effect, U is the potential and — VU is the force acting on 
the Sk, e.g., those from the boundary, gradient of the magnetic field, and the impurities. 
As for the conduction electrons, they are also subject to the effective magnetic flux 
corresponding to (j )o = 2tt h/e per one Sk in the limit of strong spin-charge coupling. 7 

In this paper, we explore the memory functions of Sks utilizing these physical prop¬ 
erties by solving the Landau-Lifshitz-Gilbert (LLG) equation numerically. In section 2, 
several elementary functions of Sks are revealed, and some models of memory devices 
are proposed in section 3. Section 4 is devoted to the summary and perspectives for 
future Sk memory. 

2. Elementary functions for Sk memory 

The model Hamiltonian for the chiral magnets defined on the two-dimensional square 
lattice is given by 

T~L — J ^ ^ vi r • (n r+ex T n r - \~ey) 

r 

-\- ^ ^ (Tlr X Tl r -i rex * D\ r ~~b Ti r X Ti r _|_ e ^ • D 2r ^ 

r 

^ h r n z r . (2) 

r 

Here we take the lattice constant a as the unit of length. The DM vectors JD\ and D 2 
specify the helicity (direction of vortex) and type (Bloch or Neel) of Sk, i.e., D\ r = De x 
and D 2 r — De y for Bloch Sk and D\ r = —De y and D 2 r = De x for Neel Sk, where 
e x and e y are the unit vectors in the x- and y- directions, respectively, and the sign 
of D is responsible for the helicity. The combination of the exchange interaction J and 
D produces the single-g helical state with q = tan -1 (D/J) under zero magnetic field. 
Usually DC J, and hence the helix period y = is much longer than the lattice 
constant, and hence the continuum approximation is justified. The critical magnetic 
field h c 32 ~ M separating the SkX state and the F state in the ground state is h c « 
0.78 x ( D 2 /J ). In the numerical study shown below, we apply a bias magnetic field h 0 
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being slightly larger than h c where the perfect F state is the ground state and therefore 
an Sk is defined as a meta-stable object. 

We study the creation, annihilation and drive of an Sk for various parameters in 
the Hamiltonian Eq. (2), system size L and boundary conditions, using the stochastic 
Landau-Lifshitz-Gilbert equation given by 


dn r 
d t 


= 7 


m 

dn r 


+ hr(r){t ) 


xn r +an r x 


dn r 


d t 


(jr ' V) Tlrd-fd [fX r X (j r • V) Ti r . (3) 


Here the Gaussian noise held hx(r)(t) = ( hr( r ),x(t ), hr(r),y(t), hr(r),z{t )) with the statisti¬ 
cal properties (. hr( r ),v{t )) = 0 and (hT( r )^(t)hT(r'),u(t')) = 2k B T(r)aS t ^8(r — r')8(t — t') 
(fx, v — x, y, z ) is introduced to examine the effect of local heating by a spatially de¬ 
pendent temperature T(r) on the Langevin dynamics of the system. We use the Heun 
scheme to solve this equation . 35,36 The last two terms represent the spin-transfer-torque 
(STT) induced by the (spin-polarized) electric current density j r . Below, the time and 
j = \j r \ are measured in units of 1/(7 J) and 2 eqJ/pa 2 (p: polarization of magnet), and 
the units are typically 1/(7 J) ~ 0.7 ps and 2 epJ/pa 2 ~ 1.0 x 10 13 A/m 2 , respectively, 
if we assume 7 = g s pB/fo (<]«'■ electron spin g-factor, p B : Bohr magneton), J ~ 10 -3 eV, 
p = 0.2 and a = 5 A. 


2.1 Write and Erase methods 

2.1.1 Heat-control Sk memory. 

The Sk creation/annihilation in the F background is a topological transition of the mag¬ 
netic texture associated with the singular configurations of magnetic moments. These 
singular configurations can be realized by local heating 37 (see Fig. 2 and the Movie 1 in 
the Supplementary Informations). Figure 2 (a) is a closeup (100x100 in L = 300 x 300) 
of the F state where all magnetic moments are perpendicular to the plane. A parameter 
set {J = 1.0, D = 0.15, h r = h 0 = 0.02, a = 0.01} and periodic boundary condition 
(PBC) are used. The Sk creation is achieved by applying a local heat corresponding to 
the temperature k B T = 1.5 J (~ T c ) with a duration th eating = 200 inside the red circle 
(15 in radius) in Fig. 2 (a), and Fig. 2 (b) shows the created Sk. This heating effect is 
introduced in Eq. (3) by the Gaussian noise field h T ^{t), i.e., k B T(r){t ) = 1.5J in the 
red circle in Fig. 2(a) for 0 < t < th eating and k B T(r)(t ) = 0 otherwise. By local heating, 
the magnetic moments at around the heat-spot are strongly disordered and constitute a 
nucleation center for the Sk creation. This condition depends on the heat intensity, the 
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Fig. 2. (Color online) Writing and erasing an Sk by heat. A parameter set 
{J = 1.0, D = 0.15, ho = 0.02, fcgT = 1.5, theating = 200, a = 0.01, L = 300 x 300} and the periodic 
boundary condition (PBC) are used, (a) Initial state without Sk. A 100x 100-region in the system is 
shown. Red and blue circles represent the heating spot sizes (the radii are 15 and 20, respectively, in 
this case) for Sk writing and erasing, respectively, (b) Sk writing by local heat, (c) Initial state with 
an Sk. (e) Sk erasing by local heat. 


heat-spot size and also the Gilbert damping constant a. For smaller a, the dynamics 
of magnetic moments becomes more active. The active dynamics, however, is rather 
destructive for Sk, i.e., too much heat is not suitable for Sk writing. Conversely, the 
strong heating is useful to erase the Sk. Figures 2(c)—»(d) demonstrate the Sk erasing 
by the stronger heating (see also the Movie 2 in the Supplementary Informations): An 
Sk is prepared at the initial state as shown in Fig. 2(c). The heat intensity GT = 1.5,/ 
and heating duration theating = 200 are the same with those used for the Sk writing in 
Fig 2(a)—»(b), but the larger heat-spot with radius 20 is used (the size is shown by the 
blue circle in Fig. 2(a)). This heat destroys the Sk as shown in Fig. 2(d). 

2.1.2 Cell memory: magnetic field control. 

One can write and erase an Sk by changing the magnetic field strength as a function of 
time in the presence of the system boundary as shown in Fig. 3 (see also the Movie 3 in 
the Supplementary Informations): A parameter set {J = 1.0 ,D = 0.18, ho = 0.03, a = 
0.04, L = 50 x 50} and open boundary condition (OBC) are used. A trapezoidal mag¬ 
netic field pulse is applied inside the red rectangle R = 20 x 30 in the system (see 
Fig. 3(a)), so that the magnetic field Hr = h r€ R is reduced (the minimum is h\ = 0.01) 
during 0 < / < 3000 (see Fig. 3(e)) and for r ^ R, h r = h 0 (bias field). The re¬ 
gion R involves the upper edge of the system and the Sk appears from the edge as 












Jpn. J. Appl. Phys. 


REGULAR PAPER 


(a) 


0 2000 4000 6000 8000 

t 

Fig. 3. (Color online) Writing and erasing an Sk by magnetic field, (a) Nano-sized chiral magnet 
without Sk. This is the initial state. The red rectangle represents the region R = 20 x 30 where the 
magnetic field pulse is applied. A parameter set 

{J = 1.0, D = 0.18, ho = 0.03, h\ = 0.01, /12 = 0.05, a = 0.04, L = 50 x 50} and open boundary 
condition (OBC) are used. Time evolution of the magnetic texture in Sk writing is shown by the 
snapshots at (b) t = 1200 and (c) 4000. Time evolution of the magnetic texture in Sk erasing is 
expressed by the snapshots at (c) t = 4000 and (d) 7700. (e) Time dependence of the magnetic field 
hR in the region R. 




Figs. 3(a)—>(b)—>-(c) by the magnetic field pulse. As seen in Fig. 3(b), in particular, the 
nucleation center for the Sk, the black area, penetrates from the upper edge (see also 
the Movie 3 in the Supplementary Informations for the whole process). Note that the 
intensity of the pulse (h 0 — hi = 0.02) is two orders of magnitude smaller than J. It is 
confirmed that the magnetic field pulse with this intensity can not change the magnetic 
texture at all when R is deep inside of the magnet and far from the edge of the system. 
At the edge of the system, the magnetic moments are lying parallel to the boundary 
and the spatial discontinuity reduces the topological stability, thus the Sk writing is 
achieved by a small magnetic field pulse. 

For the Sk erasing, in the present example, a positive magnetic field pulse with the 
maximum = 0.05 is applied at 5000 < t < 8000 as shown in Fig. 3(e). Here, the Sk 
becomes unstable and is pushed out from the region R. As implied from Eq. (1), the Sk 
moves along the boundary of region R due to the magnetic-field gradient. Consequently, 
the Sk turns in the dead-end at upper-right corner of the system as shown in Fig. 3(d), 
and eventually disappears from the system (see also the Movie 3 in the Supplementary 
Informations for the whole process). For the Sk erasing, the narrow dead-end designed 
in this device structure is of crucial importance, i.e., the Sk does not disappear if the 
space between the region R and the right edge of the system is as large as the Sk size. 
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2.1.3 Cell memory: electric field control. 

Here, we examine the Sk writing/erasing by the change in DM interaction with a fixed 
bias magnetic field. Figure 4 summarizes the result using the same numerical condition 
as in Sec. 2.1.2, except for the DM interaction (see also the Movie 4 in the Supplementary 
Informations): We study Neel Sk in this calculation, i.e., D i r = —D r e y and D 2r = 
D r e x . For r ^ R, we fix D r = D 0 = 0.18J. The bias magnetic field ho = 0.03 J which 
is slightly larger than h c ~ 0.78 x (DfJ J) ps 0.025J. 33,34 The expression of h c indicates 
that the increase in DM interaction has a similar effect to the decrease in magnetic field. 
For 0 < t < 3000, the DM interaction D R = D reR is increased as shown in Fig. 4(e). 
Upon this change in Dr, an Sk is created as Fig. 4(a)—»(b)— ■ »(c) and this behavior is 
similar to the Sk creation in Sec. 2.1.2. 

On the other hand, the Sk annihilation occurs without the drift motion and the 
behavior differs from the case of magnetic-field change in Sec. 2 . 1 . 2 . For 5000 < t < 
7000, the D r is decreased as shown in Fig. 4(e). This decrease in Dr with fixed h 0 leads 
to instability of the Sk and the size of Sk is proportional to Dr/ J. Therefore Sk shrinks 
in size and finally disappears. It is noted that the position of Sk does not change in 
the erasing process in sharp contrast to the case of magnetic field change. This is due 
to the fact that the electric field gradient does not produce the force acting on the Sk 



Fig. 4. (Color online) Writing and erasing an Sk by electric field, (a) Nano-sized chiral magnet 
without Sk. This is the initial state. The red rectangle represents the region R = 20 x 30 where the 
electric field pulse is applied. A parameter set 

{J = 1.0, D 0 = 0.18, Di = 0.28, D 2 = 0.08, ho = 0.03, a = 0.04, L = 50 x 50} and OBC are used. Time 
evolution of the magnetic texture in Sk writing is shown by the snapshots at (b) t = 1200 and (c) 
4000. Time evolution of the magnetic texture in Sk erasing is expressed by the snapshots at (c) 
t = 4000 and (d) 5820. (e) Time dependence of the DM interaction Dr in the region R. 
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since the energy of the Sk does not depend on D. 38 

2.1.4 Cell memory: Current control. 

Here, we examine the Sk writing/erasing by an electric current. A previous study 39 
shows that the notch structure in the constricted geometry is effective to create an 
Sk. This geometry is also important to annihilate Sk by lower electric current density. 
Figure 5 shows the result with a parameter set {J = 1.0, D = 0.18, ho = 0.0278, ksT = 
0.0, a — j3 — 0.04} and OBC (see also the Movie 5 in the Supplementary Informations). 
As shown in Fig. 5(a), in the 90 x 50 system, we design a notch structure 10 x 50 whose 
left edge is at 16 from the left-end of the system. In the system with a notch structure, 
the spatial dependent j r is determined by the current conservation law divj r = 0, 39 
and j represents the magnitude of the electric current density at the right end of the 
system (j > 0 for the current flows from left to right). For 0 < t < 600, an electric 
current pulse with a maximum \j\ = 0.04 (see Fig. 5(d)) is applied and an Sk is created. 
The time evolution is shown in Figs. 5(a)—»(b)—» (c). Here the edge of the system plays 
an essential role again: Due to the DM interaction, a substantial in-plane component of 
magnetic moment appears along the edge of chiral magnet and gives rise to a magnetic 



Fig. 5. (Color online) Writing and erasing an Sk by electric current. A parameter set 
{J = 1.0, D = 0.18, ho = 0.0278, ksT = 0.0, a = (3 = 0.04} and OBC are used. The 90 x 50 system 
has a notch 10 x 50 whose left edge is at 16 from the left-end of the system, (a) The system without 
Sk. This is the initial state. Time evolution in Sk writing is shown by the snapshots at (b) t = 300 
and (c) 1500. Time evolution in Sk erasing is shown by the snapshots at (d) t — 10350 and (e) 10800. 
(f) Time dependence of the electric current density at the right end of the system. See ref. 39 for 
similar calculations. 
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texture being compatible with Sk at the corner of the notch. From the specific corner, 
an Sk is created with a help of electric-current induced STT. Note that the consistency 
between magnetic texture along the corner and winding of Sk is essential to the creation. 
At this specific corner, the Sk is also easier to be destroyed than other part of edge of 
the system. The Sk erasing is caused by an electric current pulse opposite to the case of 
Sk writing. In the geometry of Fig. 5, j = 0.02 with a time-interval 8550 <t< 10500, is 
applied: First, as seen in Fig. 5(d), the approaching Sk to the notch structure is trapped 
at the corner where the Sk appears in the writing process (see Fig. 5(b)). Later on, the 
disappearing Sk runs along the edge of the notch and turns in the designed narrow 
dead end at the upper-left corner of the system as shown in Fig. 5(e). Finally, the Sk 
disappears from the system (see also the Movie 5 in the Supplementary Informations 
for the whole process). The shape of notch structure is of crucial importance to design 
the stability of Sk, i.e., the topological surgery by the electric current controlled Sk 
memory. 

2.2 Drive methods 

The STT by (spin-polarized) electric current drives Sk motion 34,39-44 and the current is 
much smaller than that for the motion of ferromagnetic domain walls (DWs). 45-51 This 
advantage is, however, diminished by the confining geometry; in the narrow Sk lead- 
track, the magnitude of required electric current for Sk motion becomes as large as that 
in the case of DWs when the electric current j is applied parallel to the lead-track. This 
is due to the effect of the Sk confining potential (expressed by Vf7 in Eq. (1)) at the edge. 
In the confining geometry, on the other hand, an efficient Sk driving force emerges at 
the edge, i.e., the Sk gliding -motion along the edge by electric current in perpendicular 
direction. 42,43 Figure 6 manifestly demonstrates the efficiency of the Sk gliding-motion 
(see also the Movie 6 in the Supplementary Informations): Here, the parameter set is 
{J = 1.0 ,D = 0.18, ho = 0.0278, L = 600 x 50}. In horizontal and vertical directions, 
PBC and OBC are used, respectively. Figure 6(a) is the initial state where an Sk is in a 
narrow lead-track. The Sk moves by a uniform electric current density j = \j r \ = 0.001 
with a duration t = 30000 in parallel ((b) and (c)) and perpendicular ((d) and (e)) 
directions, respectively. In the case of parallel current, the Sk begins to move by the 
STT but the Sk confining potential in this narrow lead-track immediately suppresses the 
Sk motion, and eventually the Sk velocity becomes proportional to ~ (3/a . Therefore, 
the Sk moving-distance for (3/a = 1 shown in Fig. 6(b) is larger than that for (3/a = 0 
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Fig. 6. (Color online) Parallel and edge gliding motions of Sk by STT. A parameter set 
{ J = 1.0, D = 0.18, ho = 0.0278, L = 600 x 50} is used. In horizontal and vertical directions, PBC 
and OBC are applied, respectively, (a) Initial state, (b) The moved Sk at t = 30000 in parallel 
configuration with a uniform electric current density j = \j r \ = 0.001 and a = f3 = 0.04. (c) The same 
as (b) but {a = 0.04, /3 = 0.0}. (d) and (e) The edge gliding motion: The moved Sk at t = 30000 in 
perpendicular configuration with j = 0.001 for (d) a = /3 = 0.04 and (e) {a = 0.04, /3 = 0.0}. See 
ref. 43 for similar calculations. 


shown in Fig. 6(c). In the case of perpendicular current, however, /? is totally irrelevant. 
Figures 6(d) and (e) shows the results of the Sk gliding motion for /?/a — 1 and 0, 
respectively. Note that the Sk moving-distance in the gliding Sk motion is much larger 
than that in Fig. 6(b). The velocity of the gliding Sk motion is enhanced by a factor 
~l/a in comparison to that in parallel current case. In addition to the large STT in 
the gliding Sk motion, the perpendicular condition gives a significant flexibility for the 
Sk-device design discussed below. 

2.3 Read methods 

Now the methods to read out the presence or absence of the Sks are discussed. One way is 
to use the tunnelling magnetoresistance (TMR). Consider the tunnel junction composed 
of the ferromagnetic (F) and the skyrmionic thin films. The tunnelling conductance G 
of the junction is given by the expression 

G cc I <ka\H\n> \ 2 5(E n - E F )S(e ka - E F ) (4) 

n,k 

where E F is the Fermi energy, | n > and E n are the eigenstate and its energy eigenvalue 
of the skyrmionic system, while £ k(J and \ka > are those for the ferromagnet. The spin 
configuration in skyrmionic system affects | n > and hence G. Suppose the F-system 
is perfectly spin polarized perpendicular to the film, i.e., a = +1 in Eq. (4). In the 
absence of Sks, the skyrmionic system is also perfectly spin polarized, and G has the 
maximum value Go in this case. When N Sks are introduced into skyrmionic system, 
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the fraction = Nn^ 2 /A (£: Sk size, A: sectional area of the F and skyrmionic thin- 
film-junction) of the spins are reversed to down direction, which cannot contribute to 
the conductance since it has no overlap integral with the wavefunction in F-system. 

Therefore, the reduction of the conductance AG is estimated as 

AG 

— - Ni^e/A. (5) 

(jro 

When the area A is comparable to the size of one Sk, this ratio can be a reasonable 
fraction of unity, and hence can be a sensitive probe of the presence or absence of the 
Sks. 

The other method for the detection is to utilize the topological Hall effect 
(THE ). 52-54 As described above, the Sk acts as the effective magnetic flux on the conduc¬ 
tion electrons coupled to it. Each Sk is associated with the unit flux b 0 = '2nh/e in the 
strong limit of spin-charge coupling . 7 This effective magnetic field scatters the incoming 
wave asymmetrically between the right and left directions. When the wavenumber k 
satisfies the condition k£ = tt/2, the skew scattering angle has a peak value of the order 
of 60 degree and decreases as ~ 1 /(k£) for k£ » l . 54 This skew scattering produces 
the Hall voltage perpendicular to the direction of the current. When the periodic array 
of the Sks, i.e., SkX, is formed, an almost uniform effective magnetic field < b z > is 
present as 


<b z >= 


\Z30o 

2 A 2 


( 6 ) 


where a s = 2A/\/3 is the lattice constant of the triangular SkX. This produces THE 52 54 
replacing the external magnetic field by this effective magnetic field. 


3. Some models for Sk memory 

3.1 Cell memory 

In Sec. 2 . 1 , we discussed the Sk writing/erasing with the simulation results Figs. 2 , 3, 
4 and 5. These results also provide a set of design principles for the Sk cell memories; 
the size of the system and local area for external forces are scaled by the Sk size £, the 
intensity of the external forces are scaled by J and operation time is scaled by 1/(7J). 
If we assume J ~ 10 -3 eV, 1 /( 7 J) ~ 0.7 ps. Therefore, the Sk writing/erasing time 
At is of the order of nano or pico seconds. For example, in the electric current control 
device discussed in Fig. 5, the Sk writing energy is expressed to be AT x pj 2 x Ad 
where p is electrical resistivity, A and d are the area and thickness of the film sample. 
If we assume the parameter set , 34,39-41 {At = 200ps, p = lOOpflcm, j = lx 10 7 A/cm 2 , 
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A = (La 2 =)45nmx25nm, d = 10nm}, E cost = 0.0225fJ is estimated for an Sk writing. 
The F exchange coupling J should be increased for the room temperature Sk-hosting 
magnet. The increase in J causes an increase in E cost oc J (assuming the same p) 
because At oc J -1 and j oc J. Nevertheless, E cost is still several orders of magnitude 
smaller than 0.09pJ which is estimated for writing electric power in the best condition 
of the currently developing STT-MRAM. 55 Moreover, the Sk memory designed in the 
Sec. 2.1.3 works by an external electric field, giving rise to a principle for a further 
low-energy-cost Sk memory device. 


3.2 Slide-switch memory 

A possible device design of the Sk memory without Sk writing/erasing is shown in 
Fig. 7, where the red rectangles represent the stable areas for the Sk and the stay of the 
Sk in left or right red rectangle corresponds to ‘1’ or ‘O’ states (see also the Movie 7 in 
the Supplementary Informations). A parameter set {J = 1.0 ,D = 0.18, h 0 = 0.03, cr = 
0.04, /3 = 0.0, L = 90 x 50} and OBC are used. The local minima of the potential for 
the Sk is prepared by the spatial dependence of h r , i.e., the static (bias) magnetic field 
inside the red rectangles (30 x 30 each) Hr — 0.029 is weaker than that in other region 
/io = 0.03. For the switching between ‘F to ‘O’, the current driven motion of Sk, in 
particular the gliding Sk motion discussed in Sec. 2.2, is efficient. In the case shown 



Fig. 7. (Color online) Sliding-switch memory. The two red rectangles, 30 x 30 each in the system 
with the size L = 90 x 50, specifies the spatial dependence in static (bias) magnetic field h r , i.e., 
h r = Hr inside and h r = ho outside. A parameter set 

{J = 1.0, D = 0.18, ho = 0.03, hR = 0.029, a = 0.04, f3 = 0.0} and OBC are used, (a) Initial state, (b) 
The result of gliding Sk motion at t = 3000 by an electric current pulse with the maximum current 
density j = \ j r \ = 0.001. (c) The magnetic texture at t = 8000. (d) Time dependence of j. For 
t > 3000, j = 0. 
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in Fig. 7, the electric current pulse in y direction with the maximum current density 
j = \j r \ = 0.001 is enough to switch the state, whereas this intensity is too small to 
do it in the case that j r || e x . Note that j = 0 for 3000 < t, i.e., the electric current 
pulse is turned off before the Sk goes into the right stable area and inertia 43 in the 
Sk gliding-motion finalizes the switching (see also the Movie 7 in the Supplementary 
Informations for the whole process). 

3.3 Racing-circuit memory 

The gliding Sk motion also enables the flexible device design; the dosed circuit can be 
the lead-track for the Sk motion which would be impossible with the parallel current 
direction discussed in Sec. 2 . 2 . Figure 8 shows an example where the chiral magnet with 
the width of 50 forms a closed lead-track for the moving Sks in the 420 x 120 -rectangle- 
area. A parameter set {J = 1 . 0,11 = 0.18, ho = 0.03, a = 0.04, /3 = 0 . 0 } and OBC are 
used. An electric current j r with a magnitude j = \j r \ = 0.001 is flowing from outside 
to inside of the Sk lead-track and causes the motion of the Sk along the circuit wall. 
Here the magnetic field control is used for the Sk creation and annihilation (see also the 
Movie 8 in the Supplementary Informations), while all the writing/erasing methods of 
the Sk discussed in Sec. 2.1 are available. In the red rectangle region (20 x 30) shown 
in Fig. 8 (a), the magnetic field pulse hn(t) is applied. For an Sk-writing magnetic field 
pulse, the condition in Fig. 3, i.e., a trapezoidal magnetic field pulse with an intensity 
ho~hi = 0.02 at a duration At = 3000 is employed. By the four pulses at 0 < t < 33000, 
four Sks are created and keep moving in clockwise direction (see Figs. 8 (b) and (f)). For 
the erasing of Sk, a weak magnetic field pulse, which is applied at 81900 < t < 84900 
and whose intensity h 0 — h 2 = 0.01 is half of the Sk writing pulse in this case, is 
used. Figures 8 (c)—^(d) show the snapshots of the Sk erasing by the weak magnetic 
field pulse and Fig. 8 (e) is a closeup of Fig. 8 (c), where the Sk disappears through the 
upper edge (see the Movie 8 in the Supplementary Informations for the whole process). 
The intensity of the erasing pulses is not strong enough to create the Sk, however, 
the local magnetic field reduction brings about an attractive force to the Sk. Because 
the Sk motion follows the Eq. ( 1 ), the Sk moves perpendicular to the attractive force. 
Therefore, the Sk approaching to the red rectangle Sk jumps out of the system through 
the upper edge by the weak magnetic field pulse, as shown in Figs. 8 (c) and (e). 

The confining potential for the Sk at the edge is finite, and hence there exists a 
critical value in the electric current density j = \j\ for the gliding Sk motion. This is 
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Fig. 8. (Color online) Racing-circuit memory. A parameter set 

{J = 1.0, D = 0.18, ho = 0.03, hi = 0.01, h<z = 0.02, a = 0.04, f3 = 0.0} and OBC are used. In the 
420 x 120-rectangle-area, the chiral magnet with width of 50 forms a closed lead-track for the moving 
Sks. (a) Initial state without Sk. Electric current j r with a magnitude j = \j r \ = 0.001 is flowing 
from outside to inside of the lead-track. In the red rectangle (20 x 30), writing/erasing magnetic field 
pulse hR{t) is applied, (b) Four Sks created by magnetic field pulses are moving in clockwise 
direction along outer-edge of the circuit. For the writing of an Sk, a magnetic field pulse as in 
Sec. 2.1.2 i.e., a trapezoidal magnetic field pulse with an intensity ho — hi = 0.02 at a duration 
At = 3000 is used. Time evolution of the Sk erasing by an magnetic field pulse is shown by the 
snapshots at (c) t = 83200 and (d) 85000. The difference between the Sk writing and erasing 
magnetic field pulse is the intensity, i.e., an intensity ho — h^ = 0.01 at 81900 <t< 84900 is used for 
the Sk erasing, (e) A close-up of (c). Time dependence of hR. (g) Reset procedure: snapshot at 
t = 88300. (h) Time dependence of j in reset procedure. 


useful for reset procedure. Figures 8(g) and (h) show a snapshot of the disappearing Sks 
and the time-dependence of the electric current pulse for the reset procedure, respec¬ 
tively (see also the Movie 9 in the Supplementary Informations): For 0 < t < 80000, 
the same numerical condition as that in the previous case, Figs. 8(a-f), is used. All the 
Sks are pushed out by the electric current pulse at 80000 < t < 9000 with a maximum 
j = 0.002 and disappears through the outer edge of the system. 

4. Summary and perspectives 

Here it is worthwhile to compare the characteristics of skyrmion memory devices with 
other magnetic ones such as bubble 56,57 and racetrack 48 memories that also employ 
memory bits or domains in magnetic thin films. Table I summarizes the comparison 
of these three magnetic memories. The bubble memory 56,57 had been developed since 
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Table 1. Comparison of memory devices. 


Bubble 

Race-track 

Sk 

Size 

~ ji m 

10 ~ lOOnm 

3 ~ lOOnm 

Write/Erase 

magnetic field 

Current (Spin transfer 

Current, Heat, Electric 



trque) 

field 

Read 

Magneto-resistive sensor, 

Magneto-tunneling junc¬ 

Topological Hall effect, 


Hall sensor 

tion (MTJ) 

MTJ, Hall sensor 

Drive 

Magnetic field 

Current: 10 7 A/cm 2 

Current: 10 2 A/cm 2 


1960’s and were at market around 1980. The bubbles were manipulated in insulating 
iron garnet films by means of an external magnetic field and were detected by the 
magnetic field coming out of the bit by magnetoresistance or Hall devices. Compared 
with bubbles, the size of magnetic bits for racetrack and Sk memories are quite small 
that are necessary condition for the pursuit of high density memory. It should be noted 
that some of the knowledge obtained in the research of magnetic bubble device 56,57 are 
referred in the course of skyrmion memory research. 

When the latter two are compared, a most distinct feature of skyrmion device is 
the very low current density needed to drive magnetic bits, orders of magnitude lower 
than that required for domain wall motion in racetrack memory. 48 This provides high 
potential for low energy-cost devices. The storage capacity of memory device is mostly 
limited by the size of the magnetic bits. Skyrmion memory having, comparable to 
or smaller size than racetrack memory, has a big advantage in this regard as well. 
Another important characteristic is the variety of tools for writing/erasing and reading 
the magnetic bits. Racetrack memory solely relies on sophisticated magnetic tunnel 
junction (MTJ). Spin polarized current injection through a tunnel barrier provides spin 
transfer torque to reverse the magnetization. This could also be employed for skyrmion 
device as discussed by Fert et al.. 58 However, this method usually consumes rather 
large energy because high density spin polarized current has to be fed through a high 
resistance MTJ. For skyrmions, as discussed in Section 2, other techniques such as 
current flow in low resistance channel, local heat application, and application of electric 
field to control DM interaction, can be used to explore novel possibility for low energy- 
cost memory devices. As for reading the bit, MTJ is an excellent method and will be 
a choice for skyrmion device as well. However, other choices such as topological Hall 
effect may explore new methods to realize high performance memory device. 
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Since magnetic structures are used, all the devices are expected to have high enough 
endurance without any fatigue or imprint that are representative limiting problem for 
ferroelectric memory where atomic motion takes place. Due to topological protection 
of skyrmion, the retention is also expected to be as excellent as other magnetic devices. 
Therefore, skyrmion has significant potential for memory devices. 

The issues discussed in this paper have been theoretically verified to raise numbers 
of challenges to experimentalists. Materials choice 7 and micro-fabrication 59,60 have been 
investigated so that we can examine basic characteristics of elementary device actions 
such as write/erase, read, and drive. As have been the case for former challenges in 
semiconductor devices and magnetic storage, interplay between theory and experiments 
as well as basic device physics and application oriented developments will be the key 
for exploring novel science/technology based on skyrmions. 
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